arXiv: 1504.08139v2 [math.AP] 22 Feb 2016 


ON MASSLESS ELECTRON LIMIT EOR A MULTISPECIES KINETIC SYSTEM 

WITH EXTERNAL MAGNETIC FIELD 

MAXIME HERDA 


Abstract. We consider a three-dimensional kinetic model for a two species plasma consisting of electrons 
and ions confined by an external nonconstant magnetic field. Then we derive a kinetic-fluid model when the 
mass ratio me!mi tends to zero. 

Each species initially obeys a Vlasov-type equation and the electrostatic coupling follows from a Poisson 
equation. In our modeling, ions are assumed non-collisional while a Fokker-Planck collision operator is taken 
into account in the electron equation. As the mass ratio tends to zero we show convergence to a new system 
where the macroscopic electron density satisfies an anisotropic drift-diffusion equation. To achieve this task, 
we overcome some specific technical issues of our model such as the strong effect of the magnetic field on 
electrons and the lack of regularity at the limit. With methods including renormalized solutions, relative 
entropy dissipation and velocity averages, we establish the rigorous derivation of the limit model. 
Keywords. Multispecies; plasma physics; magnetic field; Vlasov-Poisson; Vlasov-Poisson-Fokker-Planck; 
drift-diffusion. 

AMS SUBJECT CLASSIFICATIONS. 35Q83, 35Q84, 35B25, 82D10, 76X05 


Contents 

1. Introduction 

1.1. The physical model 

1.2. Scaling 

1.3. The mathematical model 

1.4. Main result 

2. A priori estimates 

3. Existence of solutions and uniform in e estimates 

4. Compactness of the family of solutions 

5. Taking limits in equation (11.81) 

6. Regularity of the limit 
References 


1 

2 

3 

4 
6 
I 
10 
\2 
15 
19 
22 


1. Introduction 

In many plasma physics applications, the numerical simulation of full multi-species kinetic systems of 
equations can be extremely expensive in computer time. Indeed, since the typical time, space and velocity 
scales of each species differ from several orders of magnitude, it requires a fine discretization to accurately 
approximate the different scales. We refer to [ 2 ] for a discussion on these issues in the two-species Vlasov- 
Poisson case. Therefore, part of the problem is sometimes overcome by making simplifying assumptions on 
species with negligible contribution to the whole dynamic. In this paper, we are interested in reducing a 
kinetic model by taking the limit when the mass ratio between light and heavy particles, namely electrons 
and ions, tends to 0. 

The charged gas evolves under its self-consistent electrostatic field and an external magnetic field. This 
configuration is typical of a tokamak plasma 0110] where the magnetic field is used to confine particles inside 
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the core of the device. We assume that on the time scale we consider, collisions on ions can be neglected while 
for electrons, it is entirely modeled with a Fokker-Planck operator. At the formal level, an exhaustive study 
of asymptotic mass-disparate model with more involved collision operators such as the Boltzmann or Landau 
operator can be found in the review m of Degond. When the mass ratio goes to 0, our model converges 
to a drift-diffusion equation featuring a magnetic-field dependent diffusion matrix for the electrons coupled 
with the original kinetic equation for the heavy particles. Similar parabolic equations with non-symmetric 
diffusion for plasmas can also be found in [giisiEoiiiiiiii]. 

Hereafter, we start from the physical equations and propose a detailed scaling with respect to the ions 
time scale. In the dimensionless system, for the light species equation, the leading order terms with respect 
to the mass ratio are those related to the magnetic field and the collisions. Therefore, the resulting derivation 
is in the mean time similar to a strong magnetic held limit as in the papers of Golse and Saint-Raymond 
[29l 1^ . and to a diffusive or parabolic limit as in the work of Poupaud, Soler, Masmoudi and El Ghani in 
[431 125] . While the latter papers provide us with many tools to handle our own problem, some technical issues 
in our analysis are closely related to the special features of our model. In a single-species model of charged 
particles, the other particles density is usually given either as a static regular background or as a function 
of the electric potential, while our ions are only known to obey a non-trivial kinetic equation. Because of 
this coupling, it turns out that some extra analysis is needed to recover the regularity required for our limit 
system to make sense. Besides, we have to control the strong magnetic held. This is done by looking at the 
interplay between the increasing effect of oscillations and collisions in the asymptotic regime where the mass 
ratio tends to 0. Indeed, despite the fact that collisions and magnetic forces appear at the same order of 
magnitude in original and limit equations we consider, we prove, by making the most of special cancellations, 
that with other respects collisional effects provide a form of control on magnetic contributions. While the 
rest of our analysis seems rather robust the latter control — that relies on cancellations — does not seem to 
extend readily to other forms of dissipation operators such as linear Boltzmann collision operators. 


1.1. The physical model. To avoid unnecessary technicalities, we suppose that there is only one type of 
ion in the plasma of mass rrii and charge q, while me denotes the electron mass and —q their negative charge. 
The case of a plasma containing several types of ions can be treated in the same way. 

The number of particle of type a in the phase space volume dx dv at position a; G with velocity u S 
between time t and t + dt is fa(t,x,v) dx dv dt. The index a stands for the species of particles and can be 
either i for ions or e for electrons. The ion distribution function fi evolves according to a Vlasov equation and 
the electron distribution function fe follows a Vlasov-Fokker-Planck equation. The coupling occurs through 
the Poisson equation that relates the electric held to the densities. The equations of the model, written in 
physical units, are the following 

dtfi + V ■ + ^' A Rext) ■ V„/i = 0, 

mi 

dtfe + V ■ V,/e - ^ + VA Hext) ' V„/e = • {vfe + — V,/e), 

me Ccol rUe 

- eo/Ax4> = q{ni - Ue), 

where tcoi is the characteristic time between two collisions, is the Boltzmann constant, 9 is the average 
electron temperature and Sq is the dielectric constant. The Poisson equation involves the macroscopic densities 


Ua 



Va G {i, e}. 


As mentioned before the typical time scales of ions and electrons largely differ due to the smallness of the 
mass ratio me/rui. In applications involving magnetic conhnement fusion, ions are the particles of interest and 
approximations are made on the electron distribution function fe in order to simplify the model. A common 
reduction supposes that the macroscopic electron density is given by the Maxwell-Boltzmann density 


UMsit^x) = C{t)e , 
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where C{t) is a normalization function. The derivation of the latter from (jl.ljl . with Sext = 0, is discussed 
in [3] and obtained in [S] for a one species Vlasov-Poisson-Fokker-Planck model. To our knowledge, the case 
of a magnetized plasma has never been treated before. We stress that the presence of a strong magnetic field 
modifies even formal computations. 

The goal of the present paper is to derive a fluid equation for electrons, when the mass ratio tends to 0. 
The evolution of their dynamic shall obey an equation on the macroscopic density. In the next paragraph 
we write the system in a consistent dimensionless form to receive the fluid model in the asymptotic regime 
of massless electrons. 


1.2. Scaling. We denote by L the characteristic length of the system, to the characteristic time and Va 
the thermal velocity for the species a G {i,e}. For any other physical quantity G, we denote by G the 
characteristic value of G and G' the dimensionless quantity associated to G so that G = GG'. We assume 
that the plasma is globally neutral, which means that 

hi = he =: N, 

and that the characteristic temperature (or kinetic energy) of each species are equal. A plasma satisfying 
the latter hypothesis is called a hot plasma [1] and satisfies with our notation 

TOiP/ = rUeV^ = ksB. 


The new unknowns of the system are then defined by the following relations 

, , N , f t X V \ , . T,T if t x\ 

= i?/" (h- 1’ j ’ I j • 

iliiUx) - (1, ^). 

The dimensional analysis of HH) introduces several important physical constants of the system, namely, 
for each species a G {f, e} 


Ad 


IsokBO 

q^N ’ 


fi») 

'-p 


hR fic^) 
V ’ ^ 


rrig 

q^eKt 





which are respectively the Debye length, the plasma time, the cyclotron time and the Larmor radius. More 
details on these constants can be found in the physics literature (see 12 SO]). We mention that the first two 
are typical scales of the electrostatic effects while the last two are related to magnetic phenomena. Since the 
goal is to perform a model reduction for the electron dynamic, we choose a scaling relative to the typical ion 
time scale. In particular, it means that we choose 



From some of the original physical constants and characteristic quantities arise dimensionless parameters of 
the system, namely 


5 = 


Ad 


V = 


ksO' 


= 


e = 



AM.Bext 

In a tokamak plasma these parameters would all be small. Dealing with the asymptotic (5 —?> 0 is called the 
quasineutral limit and has been first investigated for the Vlasov-Poisson system by Brenier and Grenier in 
[niisiiin] and developed recently by Han-Kwan, Hauray and Rousset in [32132131 ] • The second parameter rj 
is called the coupling parameter, for it measures the importance of the electrostatic effects with respect to the 
thermal agitation. The third parameter fi compares the Coulomb (electric) and the Laplace (magnetic) forces. 
The limit ^ 0 corresponds to the case of a strong magnetic field and has been studied for single-species 

plasma by Frenod and Sonnendrucker in [26] and Golse and Saint-Raymond in [22 02 ■ This asymptotic 
is called the gyrokinetic or drift-kinetic approximation in the mathematics literature. The last parameter 
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e, quantifying the mass ratio, is the main concern of this paper. We are interested in the limit of massless 
electrons, namely e —> 0. The dimensionless equations write 


( 1 . 2 ) 


' dtf[ + V ■ V,/' - • V„/' + \v A Bext) • V,/' = 0, 

^ \ ■ V,/' + • V„/' Bext) • V,/' = • ivf' + V„/') 

e £ /i£^ tcol 

^ - r]S‘^A^(l>^ = n[ - n'g. 


The link between the different time scales and dimensionless parameters follows from the previous relations 
and writes 




= = ^^o, 


f, = hi = 

e r] 


Finally we consider the regime where 

to _ 1 

tcol £^' 

which means that the ratio between collision time and the observation time is of the same order than the 
mass ratio. This can be derived considering that the Knudsen number, namely the ratio between the mean 
free path of the particles and the observation length L is of order e. More details on this part of the scaling 
can be found in the paper of the physicists Petit and Darrozes [42j and in the work of Degond and Lucquin 
[HIIIIIS] and more recently Graille, Magin, Massot and Giovangigli [HI 113 HE]- This makes the collision 
the leading order term with the magnetic field in the electron equation. In its analysis, the resulting scaling 
slightly differs from the usual parabolic or hydrodynamic scaling of Vlasov-Poisson-Fokker-Planck since the 
magnetic field introduces fast oscillations when e tends to zero and it is a priori not clear that the latter can 
be controlled by the dissipation effect due to collisions. 

By taking weaker collisions, say to/^coi = one formally derives a Maxwell-Boltzmann density in the 
direction parallel to the magnetic field and a guiding center model in the perpendicular plane when e —> 0. 
This limit model features interesting and well-known physical phenomena. Its derivation and analysis will 
be investigated in future work. 


1.3. The mathematical model. From now on, we keep fixed and thus, for writing convenience, we 

assume that S = rj = n = 1. The rescaled two species (Vlasov)-(Vlasov-Fokker-Planck)-(Poisson) kinetic 
system with external magnetic field writes 

f dtf! + V ■ • V,/f + {vAB)- = 0, 

(1-3) < edtf! + V ■ - l(u A B) • V„/| = 

[ - A^<jf =nl - nl- 

where we have dropped primes and marked the dependency on e with an exponent. The Fokker-Planck 
operator on the right-hand side of the second equation of (lOJ is defined by 

Cl : + 

The external magnetic field is of the form 

B(t,x) = {bi 62 bs)^ , 

in the canonical euclidean basis of K^. The magnetic field B is assumed to belong to L°°{0, T; L“) and to 
satisfy the Maxwell-Thompson equation Va, ■ B = 0. Its effect on the particle is given by the Lorentz force, 
appearing in the kinetic equations as ’’(uAB) • Vi,”. The Gauchy problem is completed with initial conditions 
on distribution functions 

f^(t = 0,x,v) = VaG{i,e}. 
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The Poisson equation can be reformulated using its fundamental solution $, in dimension 3 


(1.4) 


$( 0 :) 


1 

47r|a:| 


The potential is then given by 


(1.5) 


= 4 > *3. rig 


$(a; - y)n’^{t,y)dy, 


where n® denotes the total charge density and is given by 

7f(t,x) = nl{t,x) — nl{t,x). 


For later use, we also introduce the current density 

fit, x) = x) - jl{t, x), where ji = [ vf-dv and je = - 

v/r 3 e 


vfldv. 


and the rescaled uniform Maxwellian, in dimension 3 


( 1 . 6 ) 


M{v) = 




(27r)3/2 


The existence of global weak solutions for the Vlasov-Poisson system, namely the first and the last equation 
of (ll.3|l with a given electron background rig, was first established by Arsenev in [T]. DiPerna and Lions 
improved conditions on the initial data for this result in [22j and adapted their renormalization techniques for 
the Vlasov-Poisson system m- Concerning classical solutions, global existence in dimension 2 was proved 
by Ukai and Okabe for well-localized initial data in [IHl- The case of dimension 3 was considered by Bardos 
and Degond in [3] . 

The theory for the Vlasov-Poisson-Fokker-Planck system, consisting of last two equations of (131) with 
a given ion background nf, is also well-known. The existence of global weak solutions is proved in [171 US 
and results on the existence and uniqueness of smooth solutions were first obtained by Degond in m and 
generalized by Bouchut in [7]. Here, we shall consider DiPerna-Lions renormalized solutions of the Vlasov- 
Fokker-Planck equation as introduced in [^. Our choice is motivated by two observations. First, the 
presence of the friction part of the Fokker-Planck operator precludes any hope for a uniform (in time and in 
e) LP „ bound of /| in the small e limit when p > 1. As Poupaud and Soler proved in m, a uniform in s 
bound can be obtained for the adequate weighted norm but only on a finite time decreasing with p, which 
is not quite satisfying. On the other hand, knowing merely that /| lies in ^ and that belongs to 

is insufficient to give even a distributional sense to the product fg'^x4’^ appearing in the electron equation of 
C3D. The concept of renormalized solutions proposes to replace the direct consideration of the troublesome 
equation for /| by a family of meaningful equations for j3{f^), fd ranging through a sufficiently large class of 
smooth functions. 

For the coupled system dni), we consider weak solutions of the ion equation and renormalized solutions 
of the electron equation, the electric potential being given by (II3D. 


Definition 1.1. We say that a triplet (/f,/|, (/''^) is a solution on [0,T) to the Cauchy problem ()1.3I) with 
initial data f™, /™ if it satisfies 

1) ff G L°°(0,T;Li /f G L“(0,r;Li „); > 0 almost everywhere for a G {i,e}, 

and, for every G 'D{[0, T) x R®) 

3) The mapping t >->■ //gg ipffdvdx is continuous and the first equation of (ll.3|l holds in the sense of 
distributions on [0,T) x R® with the initial condition f™. 

4) For every function /3 G C^(R+) satisfying 

\P{u)\<C{^+l), \^P'{u)\<C, \up"{u)\<C, 


5 




for some C > 0, the mapping t i—>■ (pl3{f^)dvdx is continuous and /3(/f) satisfies in the sense of 

distributions on [0, T) x R® 

(1.7) edtPifl) + V ■ VM) + 

with the initial condition /3(/g"). 

5) For any A > 0, = \//| + AM is such that the mapping t i—>■ ffj^e (pOg^xdvdx is continuous and 

9e^x satisfies in the sense of distributions on [0, T) x K® 


( 1 . 8 ) 


£dt9e,X + V ■ Va;0E X + Vi 


V„0e,A = 


2e6i, 


e,A 




26», 


£,X 


where the Fokker-Planck operator and the term related to the magnetic field are gathered within the 
operator Ca which may be written in the following form 


(1.9) CA{f)=yv{A{t,x)vf + VJ) 

where the matrix A(t, x) is given by 

/ 1 63 -& 2 \ 

(1.10) Av = V + V A B = I —63 1 61 u. 

V &2 -bl 1 J 


Remark 1.2. One readily checks that weak formulations of 3), 4) and 5) are consistent with estimates in 1) 
and 2). We stress that we need to use two types of renormalization for the electron Vlasov-Fokker-Planck 
equation. Actually, points 1) to 4) are sufficient to define a self-consistent notion of solution for which we 
can prove an existence result. The introduction of the additional equation (na) is inspired by [2S]. It comes 
from a renormalization of the equation satisfied by fl/M with the function s 1— >■ a/s -I- A. While equation 
(EZl) provides us with the required alternative meaning for the electron equation in (11.31) . we shall pass to 
the limit e —0 in dEHl). 


I . 4 . Main result. Our main goal is to prove the convergence of solutions to m -in the sense of Definition 

II. lb towards weak solutions of the following coupled kinetic-fluid system 

' dtfi + V ■ S/xfi - ■ ^vfi + {vAB)- \7yfi = 0, 

4 ” V3, * Je — 0 ; 
je — D(Va,7Ze 

^e: 


( 1 . 11 ) 


/*( 0 , •, •) = /r and ne( 0 , ■,■) = J fTdv. 
where the diffusion matrix is given by 


D{t, x) = A ^(t, x) = 


1 


1 + lSl 


l + h\ -63 -I- 6162 62 + bibs 

bs + 6162 1 + -bl + & 2^3 


& 2 +^ 1^3 &1 + & 2^3 


l + bl 


Since B is essentially bounded, it yields 

Dg L“(0,T;L“). 

Moreover by denoting Is the identity matrix in dimension 3, one sees that 

, , D + D^ 1 

Re(D) = 


is uniformly positive definite. 
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Remark 1.3. Let us denote by £’ = —Vxi’ the electric field. In the strong magnetic field limit, namely 
\B\ —> oo, the particles follow the guiding center dynamic (see for example sminiii] in the physics literature). 
They are mainly advected along the magnetic field and transport in the perpendicular plane, called the electric 
drift, occurs at the next order in l/|il|. The corresponding parallel and drift velocities are respectively 


{E- B)B EAB 

V|| =-, Vdrift = 


\B\ 


\B\ 


The guiding center dynamic can be recovered from the limit electron equation in (11.111) . Indeed, note that 
applying the diffusion matrix to the electric field gives 

1 „ „ 1 


DE = E- 


( 1 . 12 ) 


l + l^l 


:E AB + 


l + \B\ 


:{E AB)AB, 


1 


1+ Bl 


[{E ■ B)B - E AB + E]. 


Therefore the limit equation for the electrons formally rewrites 


dtUe + Va; • [(V|| + Vdrift + O {l/\B\^) )ne - DV xTle] = 0. 


Definition 1.4. A triplet (/i,ne,^) is called a weak solution on [0, T) of the Cauchy problem (11.111) if it 
satisfies 

1) fi G L°°(0, T; „ (d L“„); /i > 0 almost everywhere, 

2) Ue G L°°{0,T-, Ll.)] Ue > 0 almost everywhere, 

3) <^ = $ *, (/ fdv - He) G L°°{0,T-Hf) and neV,(/) G Ll,{[0,T) x 
and for every ip G T>([0, T) x R®), if G T>([0, T) x 

4) The mappings t i—>■ ipfidvdx and 1 1 —>■ ifriedx are continuous and l|l.ll|l holds in the sense of 

distributions. 


Let us state the main result of this paper. We make the following assumptions on the initial data 
(Al) Va G {i, e}, /™ > 0 almost everywhere. 


(A2) 

(A3) 

(A4) 


E 

oc^{i,e} 

fi'' £ B 


fl’^dvdx = 


fL^i\x\ + \v\^ + \\na)dvdx + 

fin ^ r 1 
Je “ 1 

f^dvdx. 


' dx < +00 


where ^ *x (J fi’^dv — f fffdv). 


Remark 1.5. The case of initial data depending on e may be treated in the same way as soon as assumptions 
(EU to (IA4I) are satisfied uniformly in e. In this case one may pick as limit initial condition for (11.111) any 
accumulation point of the sequence of initial data. 


Theorem 1.6. Under assumptions (IA1|) to (IA4I) . there exists a solution {ff, fl,4U) of (11.31) in the sense 
of Definition 11.11 Moreover, for any such solution, one has, when e —>■ 0 and up to the extraction of a 
subsequence 

ff —> fi weakly-* in L°“([0,T) x R^ x R^), 
fl—> UeM strongly in L^{0,T-, 

—> (j) strongly in L^(0, T ; W^’^) for I < p <2, 


and the limit (fi,ne,4>) is a weak solution of (11.111) . 
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The parabolic limit for the Vlasov-Poisson-Fokker-Planck system with a given ion background has been 
studied in dimension 2 by Poupaud-Soler in [43j and Goudon in [30j . Masmoudi and El Ghani generalized 
these results in any dimension in [25j . using DiPerna-Lions renormalized solutions and averaging lemmas. 
The procedure they follow was introduced by Masmoudi and Tayeb to study the diffusion limit of a semi¬ 
conductor Boltzmann-Poisson system in m- Let us point out, here and later, the importance of the latter 
paper, which provides efficient tools to derive a global in time result for our own problem. A recent paper m 
of Wu, Lin and Liu treats with this method the case of a multispecies model where several Vlasov-Fokker- 
Planck equations are coupled by a Poisson equation on a bounded domain. Let us also mention |24j . where the 
author deals with the case of a self-consistent magnetic field in the context of a Vlasov-Maxwell-Fokker-Planck 
system. 

Here, we are considering a multispecies model with an external magnetic field. This brings some new 
technical difficulties in the analysis. First, the coupling with a non-collisional kinetic equation rather than a 
fixed background of charged particles makes it harder to recover regularity sufficient to give a distributional 
sense to the limit problem. Indeed, the regularity given by Vlasov-Poisson in dimension three for the ions 
impacts that of the electrons (see Lemma [6.11 and 15^ . Besides, by considering a magnetic held of the same 
order of magnitude than the collisions, it is not clear at all that one can control the fast Larmor oscillations 
at the limit especially in the three dimensional setting. As an example of this difficulty in the non-collisional 
case of Vlasov-Poisson, we refer to the paper of Saint-Raymond [33] where only the two-dimensional case 
can be treated since the dynamic along the magnetic held lines is too fast to be captured at the limit in this 
scaling. Here, we are able to control the whole dynamic thanks to the dissipative effect of the collisions and 
the orthogonality between the Lorentz force and the velocity held (see Lemma HOI and Proposition |33|)- An 
other proof of the nontrivial effect of the B held can be seen with the diffusion matrix of the limit drift- 
diffusion equation which is anisotropic, for it contains the magnetic held effects. Let us mention that such a 
model was derived in a linear setting [i.e. with an external electric held) by Ben Abdallah and El Hajj in [5| 
with a linear Boltzmann collision kernel. 

In the rest of this paper we detail the proof of Theorem 11.61 The outline is as follows. In section |21 we 
introduce natural estimates associated with (1131). These estimates will be crucial to prove our derivation 
and the existence of solution which is briehy discussed in Section [Sj After proving the required compactness 
of the family of solutions in Section |31 we will take limits in equations in Section |5| Eventually, in Section jS] 
we shall use the algebraic structure of the limit system to gain some regularity which will allow us to recover 
(11.111) in a distributional-rather than renormalized-sense. 

2. A PRIORI ESTIMATES 

The study of the asymptotic e — >■ 0 requires estimates that are uniform with respect to e. For our coupled 
system, the only natural identities providing such bounds are mass estimates (see Lemma E21), free energy 
and entropy inequalities (see Proposition l2.1l) . Let us introduce kinetic energies associated with each species 

K^{t) = Kl{t) -\- Kl(t) where Ka = ^ [[ \v\‘^fadvdx Va G {i,e}. 

^ J JRO 

The characteristic energy due to electrostatic effects is called electric energy and reads 

E'^(i) = ^ [ |Va;</>®pfix = / (j)^n^dx, 

where the last equality stems from the Poisson equation. Let us also define the entropy of each species 

SUt) = [[ fa\n{f^)dvdx Va G {i,e}. 

J Jro 

The natural energy associated with Vlasov-Fokker-Planck type equations is called the free energy and writes, 
for our system 


U^{t)=E%t)+K^t) + Sl{t). 




We also introduce the free energy dissipation given by the following non-negative quantity 

2 


1 


1 




e|2 


dvdx = 


f 

V„a/ — 
M 


Mdvdx 


Proposition 2.1 (Free energy and entropy estimates). Suppose that (/®, /|, cjf) is a smooth localized solution 
of (II.dll . One has the following ’’entropy” estimates, for all e > 0 and t € [0,T). 

• The free energy satisfies 

U^{t)+ f D‘’’{s)ds = U{0). 

Jo 

• The ion entropy satisfies 


sfit) = s!{0) 


Proof. Multiplying the first two equations of dOD by |np/2 and |up/(2e) respectively, integrating in v,x 
and summing the two equations yields, up to an integration by parts, 

(2.1) ^ J ‘ = - ^ JJ V ■ {vff + Vyff)dvdx 

The continuity equation is obtained by integrating the fist two equations of with respect to a; and 
summing the resulting equations after dividing the electron equation by e. It reads 

(2.2) dtn^ + V, • / = 0 

Then using (12.21) and the Poisson equation, we can rewrite the second term in (12.11) 

(2.3) J ■ fdx = J \\7,y(j)^\'^dx 

Hence we get an energy estimate from (EH) and E31), namely 

(2.4) + JJy. (^/l + VyfDdvdx 

The entropy equations are obtained by multiplying the first two equations of (11.311 by In ff 1 and In f^ + 1 
respectively and integrating in x, v. It yields 

(2.5) = 0, 

( 2 . 6 ) = ^ + ^vfDdvdx, 

where we performed an integration by part of the right-hand side of the second equation. Equation ()2.5I) 
provides the ion entropy estimate. By summing (1^ divided by e and (1^ . we obtain the announced 
estimate up to an integration in time. □ 


The other type of natural a priori estimate for Vlasov-type equations is the conservation of norms. 

Lemma 2.2 {U’ norms). Suppose that iff, ff,(p^) is a smooth localized solution of (11.311 . For all e > 0 and 
tG[0,T), 

• The distribution function of electrons satisfies 

ll/e(^,-,-)llLi,„=||/:iLi.„ 

• The distribution function of ions satisfies 

ll/f(i,-,-)llLS„ = ll/rilLS„, VpG[l,oo] 
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Proof. Keeping in mind that distribution functions are positive, the integration of the Vlasov-Fokker-Planck 
equation in (ITSD with respect to x and v provides the first estimate. Let us multiply the ion equation of 
m by ifiT tv and integrate in x and v to get 


Hence H/fUiP „ 



p 


= 0 . 


is constant. Letting p go to infinity gives the limit case. 


□ 


As we are working on an unbounded domain in space, we need to control space moments of the distribution 
functions to ensure that no mass can be ’’lost” at infinity. It reduces to controlling current densities as shows 
the following estimate. 

Lemma 2.3 (First moment in space). Suppose that is a smooth localized solution of (11.311 . For 

all e > 0, t G [0, T) and a G {i, e} 

[[ \x\ff,dvdx= [ f ^■fadx+ff \x\fffdvdx 
J Jr<> Jo Jr3 \x\ J Jrb 

Proof. Multiply the first two equations of (11.31) by |a;| and integrate in x, v and t to obtain the result. 

□ 


3. Existence of solutions and uniform in e estimates 


In this section, we give an existence result for m- The a priori estimates of the previous section are 
necessary to build these solutions, by a mollification procedure. Let us mention that this result follows from 
single-species cases. Indeed the coupling between the kinetic equations of dini) is weak in the sense that, 
becanse of the form of the Poisson equation, it is possible to isolate the contribution of each species in the 
electric field The addition of the magnetic field term only cause minor and harmless modifications 

to usual proofs as it is linear and does not alter a priori estimates. For the Vlasov-Poisson part, the theory 
of Arsenev may be applied. We refer to Theorem 1.3 and 1.4] for details. The Vlasov-Poisson-Fokker- 
Planck part of (irsD may be handled with the DiPerna-Lions theory of renormalized solutions mm]. Some 
technical details may also be found in [S] and miMiiMi on bounded domains. 


Proposition 3.1. Under assumptions (EH) to (IA4I) . the system C3I) admits a solution in the sense of 
Definition \l. 11 In particular, the solution satisfies the continuity equations in the sense of distributions on 
[0, T) X K® 

(3.1) 5tn® -f Va; • = 0, VaG{i,e} 


Moreover, the following estimates hold, uniformly in t G [0,T), £ > 0 and a G {i,e}. 


(3.2) 

W{t)+ [ 

Jo 

Dts)ds < [/“, 

(3.3) 

sat) < St 


(3.4) 



(3.5) 

Vp G [l,oo]. 


(3.6) 

ff \x\pdvdx<f f ^ 
J Jrb Jo Jr^ FI 


■ fadx + 

and the distribution functions are non-negative almost everywhere. 


R6 


\x\Pffdvdx. 


From estimates obtained in Proposition [ST] we infer uniform in e estimates that will allow us to take limits 
in the following sections. Let us first give a name to the particular solutions satisfying these estimates. 
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Definition 3.2. Any triplet which is a solution of the system (jl.3|l in the sense of Definition 

o associated with initial datum satisfying (El) to (j A4j) . Rnci itsolf s&tisfyiiig cstinicLtcs of Proposition 13.11 
is called from now on a physical solution of (O. 


Proposition 3.3. A physical solution of (II.3p satisfies the following properties 

(a) Control of current densities: 

(3.7) Wjthi <C{T)+K!, 


(3.8) 


\\jl\\^,<C{T) + ^DC 


(b) Uniform bounds on the free energy and moments: 

(3.9) V [[ ff^{\x\ + Ivl"^ + \Iti ff^\)dvdx + [ \V dx + [ {s)ds < C{T) 

-^0 

(c) Consequences of (b): 


(3.10) 


\\ji\\L^(0,T-,Ll) + ||je||Li(0.T;Li) + II Vj, i/Tf 11^2(0, T;L2„) 
+ -IK?; A i?) • V„/g 11^1(0,T;Li „) < C'(r) 


for some constant ClT) independent of e and t G [0, T) 

Proof, (a) First we can control the ion current density by decomposing the velocity space in the following 


way, 


Mhi < 


[[ \v\ffdvdx+ ff \v\ffdvdx 

JJ\v\<2 JJ\v\>2 


kl>2 


< 2||/“||L=(0,T;Li J + ;^ // \v\ ffdvdx. 

^ J JR6 

We can conclude with (|A2p . 

The electron current density can be controlled by the free energy dissipation. Indeed, following the idea 
in [531 Equation 2.21], it writes 

■ ^ ‘‘{v^/f!+2y.^)y7!dv 


Je = 


1 , 


1 


1 


and thus 

Wfehi < ^ll/elUi,„ + -f\^’f!+'^vf!\'^dvdx 

and the desired result follows using (|A2p . 

(b) The key arguments here are the entropy estimates (13.2p and (|3.3p . Since distribution functions are 
non-negative, we decompose the free energy in the following form 

(3.11) + Sl+- Sl_, 

where we define, for a G {?, e}, 

fa fadvdx, 


ce _ 


= 


/a In fadvdx, 


with ln''"(s) = max{ln(s),0} and In (s) = max{—ln(s),0} for s > 0. By applying the same arguments as in 
m Lemma 2.3] and estimate (13.6L one can get the following bound on the negative part of the entropy 


( 3 . 12 ) 


St - <C+^{Kt + ||jellLi(o,t,Li)) + ^ // l^lff^dvdx, 


for some positive constant C. By inequality (j3.8|l on the electron current density and (18.1211 

^ {kI + D^{s)ds^ + CiT). 

Now with estimate (USD and the decomposition of the free energy (13.111) . one can conclude that K^, E^, 
J*D‘^{s)ds, SI _ and 5'|_|_ are uniformly bounded in e and t. Replacing the index e by i, inequality (13.1211 
holds true for the ion related quantities and the bound (13.71) on the ion current density and estimate (|3.3I) 
on the ions entropy give the boundedness of Sf _ and 

(c) The estimates on current densities follow using (a) and (b). Now, as in [521 Corollary 5.3], the following 
computation 

< J ' ^^J'edvdxdt 

< + ^ II fj^dvdx, 


provides the third bound. Finally we can control the term related to the magnetic field by noticing that, 
since {v A B) ■ v = 0, 


nAi?)-V„/|||ii(o,T;Li,j = ^ I ll\ivAB)-i2S7,^ + vyf!)yf! 


< 


/o 
\B\\l 


dxdvdt 


— / D'^dt + T sup 
£ Jo [O.T) 


We conclude with the estimates in (b). 


□ 


Remark 3.4. Actually, the estimate on the ion current density can be largely improved using a classical 
moment lemma (551 Lemma 3.1]. The latter gives that {ji)e is uniformly bounded in L°° {0,T; Lx^^). 


4. Compactness of the family of solutions 

From the estimates of the foregoing section, we can infer some compactness. When e tends to 0, estimates 
(I3.5|l and (I3.9|l give, up to the extraction of a subsequence, 

(4.1) /f —>• fi weakly-* in L°“(0, T; for p G (1, oo], 

and 

(4.2) VxJ)'' -A Vxfj) weakly-* in L°°{0,T; 

where the limit is a gradient because it is irrotational in the sense of distributions. 

Lemma 4.1. Families of physical solutions of (11.311 satisfy the following properties 

(a) {fi)^ and (/|)g are weakly relatively compact in L^([0,T) x R®). 

(b) (nf)g and (n^)^ are weakly relatively compact in L^([0,T) x R^). 

(c) fl — n%M —>• 0 in L^([0, T) x R®) when £ —>■ 0 


Proof. The first two assertions follow from the Dunford-Pettis theorem. 

(a, b) Let us define W : x e“l“l/(87r), A4 : (x,v) >->■ W{x)M(v) and : m > 0 >->■ M(ln'''M-|- 1). Note that 
W and Ad are of integral 1 and T is a convex, non-negative, increasing function satisfying 


lim 

U—¥00 


^(u) 


u 

12 


= H-(X). 











Let a G {i, e}. One sees that by the Jensen inequality 


^i^)Wdx = 


< 


< 






M 

J 06 

M 


/^In 


Mdvj Wdx 
Aidvdx 
dvdx 


J O! 

M 


+ ‘S'S - + Ka + j j \x\fadvdx + C j j fadvdx, 


for some constant C > 0. Because of the uniform estimates (iJ.sp and (ESI), we hence get by the de la Vallee 
Poussin lemma the equi-integrability of the bounded families 

{fJMh C L\0,T;LHMdvdx)) 

and 

{n^JW},cLHQ,T-,L\Wdx)). 

This yields the announced weak compactness by the Dunford-Pettis theorem. 

(c) Using the log-Sobolev inequality (see [33]), we get an upper bound for the relative entropy of the electrons 


with respect to a local Maxwellian 


/I 


'M 


dvdxds < 2 


Mdvdxds 


< C{T)e^, 

where the last inequality comes from estimate (EH). By the Cziszar-Kullback-Pinsker inequality (see [H 
Theorem 3.1 and Section 4]), we have 

nt / r r \ 2 / r \ ft 


JJ l/g — nlM\ dvdxj ds 


< 2 sup 

[O.T) 


nldx 


/I In 


f! 


IM 


dvdxds 


□ 


< C{T)e^, 

which yields the L^(0, T; convergence hence the expected result. 

The results of Lemma I3TT] are not sufficient to take limits in equations. We need to gain strong compactness 
to deal with non-linear terms. While there is no hope of doing so with the distribution functions, the 
particular averaging properties of Vlasov type equations allows us to get additional compactness results in 
space on the macroscopic densities. Indeed, by the dispersion property of the v ■ transport operator, 
one gains regularity on velocity averages of the distribution function. Besides, time compactness stems from 
the continuity equations EH) and the uniform bounds on current densities. We shall apply the following 
averaging lemma that is adapted from the famous DiPerna and Lions results in |23j . 

Lemma 4.2. Let {hL)^ be a bounded sequence in L^(0, T; ^)) satisfying in the sense of distributions 

edth^ V ■ Vxh^ = Hq + Vv ■ h\ 

where (hg)^, (hl)^ are bounded sequences in L^{0,T; Then for all tp G T>(]R^), 


J {Tyh^ — h^)'ipdv 


0 




when J/ —>■ 0 uniformly in e, where Ty is the translation of vector y in the x variable. 
Proof. We refer to EZl Appendix 2] for a proof. 


□ 
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Let us highlight that the following result is crucial to the rest of the proof of Theorem 11.61 and that in 
particular, point (b) in ProDOsition ld.dl denends upon the control of the magnetic leading term of the equation 
by the entropy dissipation. 

Proposition 4.3. Families of physical solutions of (11.311 satisfy the following strong compactness properties 

(a) is relatively compact in L^(0, T; for 1 < p < 2, 

(b) (n|)g is relatively compact in L^([0, T) x K^). 


Proof. The result of (a) stems from applying the Aubin-Lions-Simon lemma [451 Theorem 5, Corollary 4] to 
the family of electric helds and the proof can be readily adapted from m Proposition 3.3 3)]. 

(b) Let us define 

There exits > 0 such that, for all u > 0 

< min(i,M), (1 + ^)\Ps{u)\ < 1 + ^, \u/3s{u)\<^. 

In particular, one checks that Ps satisfies the requirements of a renormalization function for the Vlasov- 
Fokker-Planck equation and hence (ini) holds. Now set 

h^ = W!) 

K = -v./i • A B) ■ v.mm 

h\ = -vA^Ps{fl) + 

With the estimates of Proposition 13.31 one can show that the sequences (^•o)e and {h\)^ satisfy the 

hypotheses of LemmaWe refer to [23 Proposition 6.1] for the details. In our case, we additionally need 
to check that the magnetic field term satisfies the bound. Indeed, using that {v A B) ■ v = 0, we have 




-{vAB)-W.iPsif!)) 

{v A B)^- i(2V.v^ + v^/f!)P's{f!) 


ibO.T;Ll J 




IL“(0,T;L“) Ill'll /e' 
Then for any fixed S and any ijj G BfEP) 


II 

vf! + v„/| 

L“(0,T;Li „) 



ll/®5(/e)llL~(0.T;L“„) • 


(4.3) 


(Ty^sif!) - Ps{fl))ll)dv 


L^(0,T-,Ll.v) 


-A 0 


when 2 / —>■ 0 uniformly in e. The L°°(0,T;L^((l + \vp)dvdx) uniform bound on ff from (13.41) and 
allows us to extend this to ipiv) = 1. Now, we can also take the limit d —?> 0 uniformly in e, using the 
equi-integrability of the family (ff)e- Indeed, since, for any m > 0 

0 < Ps{u) < u and \Psiu) — m] < Su^, 

one has 

(4-4) ll/e-/35 (/|)IIli([o.t)xR6 < 2 /" [[ ffdvdxdt + 6 [ [[ \fl\^dvdxdt 

Jo JJfi>M Jo JJif;<M^ 
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for an arbitrary M > 0. The first term is 0(l/ln(|M|)) because of the uniform bound on the entropy Sg. 
Using the uniform bound on the mass, the second term is seen to be 0{5M). Take M = to conclude. 
Therefore, using (14.311 and (14.41) . one has 

when 2 / —>■ 0 uniformly in e. 

Finally, using the continuity equation em, we get a L^(0, bound on dtn%. With the uniform 

L^((l + \x\)dxdt) estimate (13.6|) on n%, this gives the relative compactness of the sequence. □ 


Using the results of Lemma 14.11 and Proposition 14.31 we get the following strong convergence results 
concerning the macroscopic density of electrons. 

Lemma 4.4. The families of physical solutions of (11.31) satisfy the following properties, up to the extraction 
of a subsequence. There exists € L^([0, T) x R.^) such that when e —>■ 0 

(a) nl —>• He in L^{0,T; L].), 

(b) y/n%^ ,/nZ in L'^{0,T-,Ll), 

(c) fl^Ug M in L^{0 ,T;Ll ,,), 

(d) = \/ ff -\- AAL —>■ y/ug -\- X\/M in L^(0, T; „), 


(f) 


;m 


Se.X 

for any A > 0. 


y/UeM 
•\/ne + A 


L^{0,T-,Ll,), 


Proof. Properties (a, b, c, d) are straightforward consequences of Lemma 14.11 fcl and Proposition 14.31 fbl. 
The last two assertions follow from (b), (d) and the bounds on y/M/dg^x and y/nffy/ug + A. □ 


5. Taking limits in equation p.8l) 

Using the previous compactness results, we can readily take limits in weak formulations of the Vlasov 
equation for ions in USD, the continuity equation for electrons (EH), and the Poisson equation in (IE3D. 


Lemma 5.1. The limits ng S L^([0,r) x R^), fi G L°“(0, T; „ ("I and S/xf) S L°°{0,T; L^) defined 
respectively in Lemma \U\ (BH and (BH satisfy in the sense of distributions, 

' dtfi + V ■ Vxfi - Vxcj) • Vyfi -\-{v AB)- Vyfi = 0, 

H” ^X ‘ je — 0; 


n„ = y fadv, Vo G {i,e}, 


for any accumulation point jg of the family {jg)e- 

At this point, we do not know much on limit^ of the electron current density j|. Our goal is now to 
characterize them. Let us introduce, 

1 


(5.1) 

and 

(5.2) 


:y/M 


(y/JI-y/^), 


Rl = - nlM) 


^From estimate 13.31 we only know a priori that, up to the extraction of a subsequence, (jDe converges to a finite Radon 
measure for the weak star topology. 
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so that the distribution function of electrons may be written as 


(5.3) 


f! = nlM+ 2erly^^M+ e^\rl\^M 
= nlM + 2eRl. 


With this in hands, one writes the electron current density as 


fe = / rlMvdv + e 


\rl\'^Mvdv 


As in [13 [H], we aim at taking the limit e —>■ 0 on the latter and characterize the limit current. We first 
gather some useful estimates on r®. 

Proposition 5.2. Let (/|)e he the electron distribution functions of a family of physical solutions of (11.31) 
and define r® by (EH). Then, the following uniform estimates hold 

(a) (rg)e is hounded in L?{Q,T-,Lf{Mdvdx)), 

(b) (e|r'gp|z;pM)e is bounded in L^{0,T; 

(c) {y/£\rl\'^\v\M)s is bounded in 

(d) (V^rg)e is bounded in Lf{{),T]L'^{Mdvdx)). 


Proof. Properties (a, b, c) can be readily adapted from [13 Proposition 5.5]. 
(d) From the definition of r® stems 


yrll'^ M dvdxdt = 


Mdvdxdt. 


The right-hand side is uniformly bounded thanks to estimate (13.91) on the entropy dissipation. 


□ 


From now on we consider a subsequence of a family of physical solutions of EH such that all con¬ 
vergence properties following from the previous compactness results hold. Let us denote by r^ the weak 
Lf{t),T-,L'^{Mdvdx)) limit of 

Lemma 5.3. When e —>■ 0, the family {fl)^ satisfies 

jl —>■ je '■= 2.^/nf / revMdv weakly in L^{0,T; L].). 


Proof. Take the limit in the following expression, coming from (15.31) . 

(5.4) = 2yJnF^ J vrlMdv + e j v\rl\^Mdv. 

The first term in the right-hand side converges weakly in L^(0,T;L;]. „). Indeed, 

^/rif strongly in L^(0,T;L^), 

and since r% —>• Cg weakly in Lf{Mdvdx)) and (t, v) v € L'^{0,T; L'^{Mdv)) 

J vrlMdv ^ J vreMdv weakly in L^(0,T;L^). 

The second term in the right-hand side of (|5.4() goes to 0 in L^{0, T ; L].) since we know from Proposition [52] 

(c) that {y/£\rl\'^\v\M)s is bounded in L^(0, T; L]J, ^). □ 
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Now we focus on the limit of equation lEHl) from Definition o to derive the last pieces of information 
we need to characterize the limit of the current density je- To do so, we introduce elements of notation 
associated with the operator £a- We denote by Lj^ the Hilbert space L'^{M.^,M~^dv) endowed with the 
scalar product 

{f,9)=[ fgM-^dv. 

Almost everywhere in t, x, 

(A(<, x)vf + V„/) = • (vf + V„/) + (u A B{t, x)) ■ V„/ 

is an unbounded operator on L\^. On 

(5.5) CaU) = • (MV„ ) + (^ A H) • V„ y/M, 

one sees that the formal adjoint of Ca in is given by 

(5.6) C*A{f) = V, • {vf + V,/) -{vAB)- V,/ = V, • {A^vf + V„/) = Cat (/). 

Proposition 5.4. The limit density and electric field —\7x4> are such that x\/^ ~ x4>\/^) G 

L^{0,T', L^). Moreover the limit electron current density je satisfies in the sense of distributions 

je — ^\/rieA (V xy/Tle x4^y/ Tie)-i 

Proof. We know from Lemma 15.31 that 

je = 2y/nl J revMdv. 

We will now take the limit in the renormalized equation dlHl) in order to get an additional equation charac¬ 
terizing re- Let us recall that (11.811 reads 


sdtOe X + V ■ A + 


= 


2ee, 


£,X 


r am 


29, 


e,\ 


By the strong convergence of 9e,\ from Lemma [4.41 Id) and the weak convergence of Vx(IA from (14.21) . the 
left-hand side of equation (11.81) converges to 


Va; • {vyj{Ue + X)M) + • (V^f/'x/(Ug -|- X)M) 

= vx/M ■ {Vxx/ne + \ - ^V^f/'x/ue + A). 

The first term of the right-hand side of the renormalized equation (11.81) may be written, using (15.21) and 
the fact that M G Ker£Ai, as 

(5.7) 

C^£,A 

Now, for any (p £ P([0,T) x E®), one has 

r = ^ JJ -^CA{Rl)Tdvdxdt := !{ + /| + J|. 
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with 




II ■■= - 


II ■■= - 


ipM 


{v A B) ■ Vi 


0e,X 

[Rl 


M 


dvdxdt. 


>0 


6 


'e,\ 


M 


{Jm) 


Ad dvdxdt, 


(^) (px/Mdvdxdt. 


where we have performed an integration by part on the Fokker-Planck part of the operator Ca defined in 
(E3) to obtain /| and /|. Mark that 


0. 


1 


'(S) 


3/2 


SO that using the definition of i?g in (15.21) . one has, for some constant C{T), 


141 < 


1 


M 


. fe 

V\ — 


dvdxdt 


n 


eix 


L°°( 0 ,r;L“^) 


‘P 


^e,X 


L“(0,T;L“„) 


< C{T)e, 


by the control on the free energy dissipation provided by estimate (13.91) . Hence /| goes to 0. To handle If 
and /f, we decompose according to, 

§ = + iKf 

Doing so, the contribution of the magnetic field becomes 

^ ^ y/UeAd 


If = 


\ 0e,X \l rZe + A 


ip(v A B) ■ Vyvldvdxdt 


[ [[ :rlx/M{v A B) ■ Vyif dvdxdt 

Jo J J 'Jrie X 


+ e 


0 J J Vne + A 

JJ A B) ■ X/yrl'/Mrl'/Mdvdxdt. 


after an integration by parts in the second term. The first and third terms of the right-hand side go to zero 
bv Lemma 14.41 ID and l5.2l lab('dl. Therefore 


4 


Jo J J '‘/Ue + X 

fT 


Jo J J \/ne + Xx/M 

The contribution of the Fokker-Planck part may be written as 

f xfn^Xd y/UeAI A 


r^x/Miv A B) ■ Vyipdvdxdt 

{v A B) ■ Vy{rf,Ad)ipdvdxdt 


I! =- 


\ 0e,X \/rie -|- A 


^ x/MX/yvl ■ \7y dvdxdt 


y/Ue + X \ 


— e 


(^MVy dvdxdt 
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As for If, the first and third terms go to 0 and therefore 

pT 

r. 


, ne 

^ lo JJ 'Jrie + XVM 

We eventually showed that, in the sense of distributions 


(5.8) 


1 




V„ • (MVi, (re)) Lpdvdxdt 

CAireM) 


2e0e,x ' y/K^TxVM 

The convergence of the second term of the right-hand side of the renormalized equation (11.81) stems from 
the strong convergence of M/9g^\ from Lemma [4.41 and the weak convergence of Vx4>^ from (14.21) . Finally, 
we receive for any A > 0 


V-/M- (^x\/ne + A - ]^yx(t>\/ne + A^ 


\/ {ne + A)Af 


CA{reM) - 


AM 


2'\/ + A)M 

By dominated convergence, one may take the limit A —>■ 0 to obtain 


V ■ V. 


vM ■ {VxV^ - = CA^TeM) 

in the sense of distributions. Since the left-hand side of the former equality is rapidly decaying in v, one 
may actually multiply the previous equation by v and integrate in the v variable to derive in the sense of 
distributions 

~ ' vCA{'i’eM)dv 


(^ xx/‘^e 2 ^^e) — 


= J C\{vM)rgdv 

= J [Ci{vM) — {v A B) ■ yv{vM)]redv 
= J —{v + {v A B))rgMdv =—A J revMdv, 

Since A(t, x) is invertible one gets the result by combining this identity with the expression of je from Lemma 

15.81 □ 


6. Regularity of the limit 

Let us summarize what we have proved so far. The triplet x4’,ne) S L°“(0, T; „ n L“„) x 
L°“(0,T;L^) X L^{0,T-,Ll.) is such that in the sense of distributions 

' dtfi + V ■ Vxfi - ^x(t> ■ ^vfi + {v AB)- Vyfi = 0, 

dt'IT'e T ^X ‘ je — 0; 

< (Vx\/Re ~ x4’x/^) G A^(o,r, L^) 

je — (V x\/’^& 2 ^ x^x/ Re) 

^x^ — Tie, 

with initial data 

/i(0, •, •) = fr and ne(0, •, •) = J ffAdv =: n“. 

Note that because of uniform bounds on {dtn%)e in and on {dtfi)e in L^{Q,T-,W~l'^), we get 

that, by the Arzela-Ascoli theorem, the limit functions t i-A f rie^pdx and t i-A f fi'tpdvdx are continuous on 
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[0, T) for any test functions ip and ijj. Thus, we do recover the above initial conditions for the limit system. 
On the other hand, usual arguments based on the convexity and lower semi-continuity of the energy and 
entropy functionals and corresponding uniform estimates prove the boundedness of the following quantities, 
uniformly in t G [ 0 , T) 


ne\lnne\dx + J \x\nedx + JJ {\v\'^ + \x\)fidvdx < C{T) 
f^\\nf,\dvdx < ll/rin/riUl., 


for some positive constant C(T). Furthermore, the following mass estimates and global neutrality result hold 
uniformly in t € [ 0 , T) 

(6-1) ll^e||Li = ||/d|Li,„ = ||/“||Li,„ = ||/riUi,„. 

(6-2) ll/illLS.„ < Wfrhljor P G (1, +oo] 

Equality is verified for the mass estimate thanks to the tightness of the distribution functions that comes 
from the control of space and velocity moments. 

The classical moment lemma m Lemma 3.1] shows that, by the boundedness of fi in L°°{0,T; L^{{1 + 
\v\'^)dvdx) n the limit macroscopic ion density has the regularity 

(6.3) mG 


We may gain some additional regularity on rig using the particular structure of the limit system. The 
procedure we set up is a generalization of Lemma 7.1 in m- Here the situation is trickier because the ion 
background rii is not regular enough to reach directly an regularity in space for Ue in order for the product 
UeVx4' to make sense. However one may first gain some regularity and conclude by a bootstrap argument in 
Lemma 16.21 


Lemma 6.1. Let Ue G L^([0, T) x be a non-negative function that satisfies 

(6.4) \/^y/nf-\-^Ey/nf=G, 

(6.5) Vx ■ E = Ui- Ue 

in the sense of distributions, where G G L^([0, T) x E G L^([0, T) x M.‘^) and Ui G 77'([0, T) x R'^). Then, 
for any p G [1,2] it holds 

Hi G LP{[0,T) X R‘^) ^UeG L^([0,T) x K'^) 


Proof. Let p G (1,2]. The first step of the proof is the renormalization of equation (16.41) . We define hereafter 
the particular renormalization function we use and which is built to recover in the end an bound on Ue- 
Let us define 7 G C°“(K+) such that 7 ( 5 ) = s on [0,1], 7 ( 5 ) = 2 for s > 3 and 0 < 7 ' < 1. Now set, for 
(5 G (0,1], 

75 ( 3 ) = 7 (<5s) Q7(<5s) + 1^ 

The derivative of the renormalization function satisfies 

( 6 . 6 ) | 75 (s)| < 1 and |s 75 (s)l < 

Equation ()6.4I) implies that x-\fnf G Ll^fi[0,T) x R^*). Let us renormalize equation ()6.4I) by multiplying it 

by 

’^x^siVnf) + 2^V^75(\AiI) = 
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One can check with (16.6|) that every term is square integrable. By taking the norm of the equation and 
expanding we obtain 


(6.7) 




+ 


JJ '^xlsiV^) ■ Ey/nl'y'^[y/^)dxdt < ||G|||2. 


We want to rewrite the third term as the scalar product of E with a gradient in order to use 
define 

( 6 . 8 ) 75(s) = [ {n's{u)fudu. 

JO 

Using (|6.8p in the third term of (1^ and dropping the first two non-negative terms yields 

JJ E ■ Vx^siy/n^)dxdt < ||G|||2. 

By using equation (16.511 , one gets after integrating by parts 


. Let us 


(6.9) 


jj Jsiy/n^) {pe - n^) < IIGIII2. 


Let us estimate ^s{y/n^)- For s G [0,1/(5], one obtains 

ls{s) = [ u{u+l)'^P~*du < [ v^^~^du= , 

Jo Jo 2p-2 


o 2 p -2 


and 


75(s) = 


(s + l)2P-2 - 1 {s+ l)2P-3 - 1 3 

-II n ^ , 

2p-2 2p-3 ^^2’ 


s — ln(s -I- 1) 


■f 3 
ifp= 2 - 


Then, one readily checks that for any p G (1, 2] there exists Gi > 0 such that for any 6 G (0,1] and s < ^ it 
holds 

75(s) > - Cl 

For s > 1/5, there exists G 2 > 0 depending only on p such that 

75 ( 5 ) < 75(3/(5) < J ^{ufP~\du = C2S^P~^ 

This provides the following estimates, for any s > 0 and G = max(Gi, G 2 , l/(2p — 2)) 

(6.10) - G^ 1 L{s<i/ 5 }(s) < 75 ( 5 ) < Gs^P"^ Vs>0. 

We now use estimate (I6.10|) in (16.9|) . First we get rid of the part involving the ion density n^. Using Young’s 
inequality, for any t] > 0 there exists some constant Cr/ > 0 such that 

75(v^)«-i < Cn\\ni\\% + r] JJ dxdt. 

< Cr,\\ni\f^p + c~p JJ 

where we used estimate (16.101) in the second inequality. On the other hand, we have 

ne^siy/n/)dxdt > ^ // nldxdt - C\\ne\\L^. 

JJ ^ J J{n^<l/S'^} 
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As a result, by taking t] sufficiently small, one gets 

jj '>T'e^{n,<l/S^dxdt < C{\\G\\L2,\\ni\\LP,\\ne\\L^), 
uniformly in (5. Taking the monotone limit <5 —>■ 0 concludes the proof. 

□ 

The regularity of rig that one may establish by the previous proof is limited by the regularity of rii. 
Nevertheless the available regularity of rii is sufficient to provide us with a termwise sense for je- 

Lemma 6.2. Limiting densities of families of physical solutions satisfy 'S/x4>y/^ S L^(0, T; ^), G 
L^{0,T; and in the sense of distributions 

je — L){^ xne 

where we recall that D = A~^. 

Proof. First by an application of Lemma l6 .1 1 with p = 5/3, we show that the source n in the Poisson equation 
—Ax4> = Hi — Ue ='■ n lies in T; fl L°°(0,T;L^) and therefore, with an L°°{0,T; L^) electric 

field, it yields Vx<?!> = Va;d) *x n. Thus, by the Hardy-Littlewood-Sobolev inequality, Vx^) G T; Lx^^'^). 

Hence, by the Holder inequality, the product Vxi/i/nl is in T; and since 

x\f^e + £ -^^([ 0 ) ^ 1 ^^)) 

this yields yjnf G T; By Sobolev embedding this gives at least .Jnf G L^(0, T; and 

since V xf> G L°°(0,T;L^), the product xf>\f^ belongs to L^(0, T; which yields the results. □ 

This completes the proof of Theorem 11.61 
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